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Abstract 

Based on the recent formulation of a general scheme to construct boundary Lax 
pairs, we develop this systematic construction for the An affine Toda field theories 
(ATFT). We work out explicitly the first two models of the hierarchy, i.e. the sine- 
Gordon (A^ ) and the A^p models. The A^p Toda theory is the first non-trivial 
example of the hierarchy that exhibits two distinct types of boundary conditions. We 
provide here novel expressions of boundary Lax pairs associated to both types of bound- 
ary conditions. 
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1 Introduction 



The Lax representation of a classical dynamical system consists in the formulation of two 
endomorphism-valued objects (in most cases matrix- or operator-valued) L and M, depend- 
ing on the dynamical variables, such that the equations of motion are contained in the 
iso-spectral evolution equation: 



The Lax matrix L therefore lives in some Lie algebra, finite-dimensional, loop algebra, or 
differential algebra depending on the specific system. It is collectively called "auxiliary 
algebra" . 

The spectrum of L, or its associated monodromy matrix if L is a differential operator, 
provides therefore a generating set of natural time-invariant candidates to be identified as 
integrable Hamiltonians. Liouville integrability of any dynamical system associated to a 
given function on this set follows from the Poisson-commutation of the elements of the set, 
guaranteed by the necessary and sufficient condition [21 [3] of the existence of a generically 
linear Poisson structure characterized by a so-called classical r-matrix, for the spectrum- 
generating operator (Lax matrix or monodromy matrix). A natural construction for M, 
given any specific Hamiltonian built as a function / on sp(L), is then available in terms of 
df, L and the r-matrix [H [2]. 

In our previous paper [I] we have examined the situation when the Poisson structure 
available for L is expressed in terms of a non-dynamical classical r-matrix plus a set of non- 
dynamical parameters encapsulated in a "boundary" matrix K obeying some purely algebraic 
quadratic equation in terms of this r-matrix. This more complicated equation (actually two 
distinct forms thereof) represents a classical version of the quantum Cherednik-Sklyanin 
reflection algebra [3 El [7] . The operator generating the Poisson commuting Hamiltonians 
then combines L and K. We have then defined a systematic construction of the M operator 
in terms of r, L and K 

Note here that the denomination of K as "boundary" matrix reflects the fact that such 
parameters encode indeed boundary effects in the Hamiltonians when the matrix L is a 
coproduct of local Z« matrices on a finite lattice % = 1, ...,N or a monodromy of local 
differential operators on a finite line. When by contrast L is a purely local Lie-algebra 
valued matrix, such parameters may be better characterized as coupling constants in a 
folding procedure (one is then considering a system on a one-site space lattice for which the 
notion of a "boundary" has no sense). 

Two fundamental situations were described in [I] : linear Poisson structure and quadratic 
Poisson structure for the generating operator. The second one is relevant to describe systems 



d t L 



L, M 
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on a lattice or a continuous line. We shall here restrict ourselves to this latter case, partic- 
ularizing it even more to degree-one differential operators L = d/dx + l(x). The generating 
operator is here the monodromy matrix of this differential operator computed between the 
two ends of the finite line for x. 

We consider here such operators for which the r-matrix is the one associated to the generic 
affine An Toda field theories. In this way we construct the associated Lax representation 
of ATFT with non-trivial integrable boundary conditions parametrized by the K matrix. 
Two "boundary" reflection equations must be considered, resp. characterized as "soliton- 
preserving" (SP) and "soliton non-preserving" (SNP) (see relevant studies at quantum and 
classical level [8] [26]). More general reflection equation may be considered given a classical 
r-matrix, through the choice of some auxiliary-space anti-automorphism, but this extension 
shall be postponed for further studies. It will be noted that this scheme automatically yields 
boundary conditions compatible with integrability. Careful evaluation of the contributions 
to the Hamiltonians and the M matrix at the edges of the x line are required to get consistent 
results. 

Our essential purpose here is two-fold. We shall first validate our general derivation by 
a comparison of the A± and Soliton-non-preserving cases with known results on low- 
dimension ATFT obtained by case-by-case analysis, basically for SNP boundary conditions 
[9]. We shall indeed show that both derivations yield exactly identical formulations for the 
boundary conditions for the dynamical fields and related Lax formulation (with a proviso, 
related to analyticity conditions, to be specified later). Similar consistency checks were also 
achieved in [I] within the vector non-linear Schrodinger context. 

We shall then provide novel expressions for boundary Hamiltonians and the associated 
Lax pairs in the yet untreated case of A*jp with SP boundary conditions, and we shall address 
some intriguing technical points that arise even in the simplest case, i.e. the sine-Gordon 
model. 

2 The general scheme 

Our analysis of classical integrable field theories with integrable boundary conditions relies 
on the study of an associated auxiliary linear problem (see [27] and references therein). Let 
us first recall the bulk case (no extra "boundary" parameter) to fix the notations and recall 
the basic structures. 

The Lax pair [28] formulation [29] of classical integrable Hamiltonian systems consists 
in defining an auxiliary linear differential problem, reading in the simplest case (order- 1 
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differential operators): 



(a t -v(x,t,\) > )v = Q 



(2.1) 
(2.2) 



U, V are in general n x n matrices with entries functions of the dynamical fields, their space 
derivatives, and possibly the complex spectral parameter A. Compatibility conditions of the 
two differential equations (12. ip . (12.21) lead to the zero curvature condition [29l 130] I3T] 
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U, V 



0. 



(2.3) 



The latter equations give rise to the corresponding classical equations of motion of the system 
under consideration. Natural conserved quantities are obtained from the monodromy matrix 



(2.4) 



T(x, y, A) = V exp{ / U(x', t, X)dx'} 

J X 

once it is assumed that U obeys the classical linear Poisson algebraic relations: 
U«(x,A), U 6 (y,/x) 



(2.5) 



(2.6) 



(2.7) 



: r ah {\ - fi), U a (x, A) +V b (y,fi) 8(x-y), 
and consequently T(x,y,X) satisfies (see [27]): 

T a (x,y,t,X), T b (x,y,t,fi)j = r ab (X-fi), T a (x,y,t, X)T b (x,y,t, fi) 
One then immediately gets: 

jtrT(A), trT(/x) | = where t(X) = trT(X) 

and thus trT(X) yields the relevant conserved quantities. 

We are now interested in implementing non-trivial integrable boundary conditions. We 
focus here on two distinct types of integrable boundary conditions: the so-called soliton pre- 
serving (SP) and the soliton non-preserving (SNP). Formulation of the two distinct types of 
boundary conditions is achieved by defining two types of Poisson structures for the modified 
monodromy matrices T. These will in fact represent the classical versions of the reflection 
algebra R, and the twisted Yangian T written in the following forms (see e.g. [6l [32], [33j [9]): 



(I) SP (reflection algebra) [5],[6] 

{7i(Ai), T 2 (A 2 )} = r 12 (Ai-A 2 )T 1 (Ai)r 2 (A 2 ) 
+ T 1 (Ai)r 21 (A 1 + A 2 )T 2 (A 2 ) 



T 1 (A 1 )T 2 (A 2 )r 21 (A 1 -A 2 ) 
T 2 (A 2 )r 12 (A 1 + A 2 )T 1 (A 1 ) 



(2.8) 
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(II) SNP ((g)-twisted Yangian) 

{^(AO, T 2 (A 2 )} = r 12 (A 1 -A 2 )T 1 (A 1 )r 2 (A 2 )-r 1 (A 1 )T 2 (A 2 )r^ 2 (A 1 -A 2 ) 

+ ^(AOr^CAi + A 2 )T 2 (A 2 ) - T 2 (A 2 )r* 2 1 (A 1 + A 2 )T 1 (A 1 ) (2.9) 

In most well known physical cases, such as the AjJ-_j r-matrices r^ 2 t2 = r 21 , hence all the 
expressions above may be written in a simpler form. 

In order to construct representations of ( 12. 81) . ( 12.91) yielding the generating function of 
Poisson-commuting Hamiltonians realizing the integrals of motion for a new classical inte- 
grable system, one now introduces non-dynamical representations (K^) of the algebra M 
(T). The non-dynamical condition: 

{jff(Ai), i^(A 2 )}=0 (2.10) 

transforms ( 12.81) . H 2 . 9 1) into algebraic equations for K ± . We consider then any bulk mon- 
odromy matrix T with Poisson structure ( 12.61) and we define in addition: 

f (A) = T~\-\) for SP, f (A) = T'(-A) for SNP. (2.11) 

One expects that a more self-contained formulation may involve a more general anti-automorphism 
of the auxiliary algebra where T lives. A corresponding reformulation of the relevant Poisson 
structures should also be proposed in that framework. We shall leave this generalization for 
later studies. 

Generalized "monodromy" matrices, realizing the corresponding algebras M, T, are finally 
given by the following expressions [6l [9] : 

T{x, y, t, A) = T{x, y, t, A) K~{\) f{x, y, t, A). (2.12) 

The generating function of the involutive quantities is defined as 

t(x, y, t, A) = tr{K + (X) T(x, y, t, A)}. (2.13) 

Indeed one shows: 

{t(s,j/,t,Ai), t(x,y,t,\ 2 )} = 0, A 1) A 2 eC. (2.14) 

The systematic Lax formulation in the case of open boundary conditions is described in 
[I]. More precisely it was shown in |1] that 

In t(X), U(x,/i)| = d V (^ A ^) + [v(x,A,/i), U(x,/x)l (2.15) 
J ox L 
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where we define: 



(I) SP 

V(x,A,/i) = r\X) tr a (K+(A)T a (0,x,A)r afe (A-/i)T a (a;,-L,A)if-(A)T ( r 1 (0,-L,-A) 
+ A+(A)T a (0, -L, A)JC(A)T a " 1 (x, —L, -A)r 6a (A + fi)T~\0, x, A)) , (2.16) 



(II) SNP 

V(x,A,/i) = r\X) tr a (^K^(X)T a (0, x, A)r o6 (A - fi)T a (x, -L, \)K~(\)T* a (0, -L, - A) 

+ iC(A)T a (0,-L,A)iC(A)T^ (2-17) 



Particular attention should be paid to the boundary points x — 0, — L. Indeed, for these 
two points one has to take into account that T(x, x, A) = T(x, x, A) = I. 

(I) SP 

V(0,t,A,/i) = r\\) tr a (K+(A)r ab (A-/i)T a (0,-L,A)JC(A)T- 1 (0,-L,-A) 
+ A+(A)T a (0, —L, A)A7(A)T a " 1 (0, — L, -A)r 6a (A + /i)) 
V(-L,t,A,/i) = r x (A) tr a (A; + (A)T a (0,-L,A)r a6 (A-/i)A-(A)T ( T 1 (0,-L,-A) 
+ A+(A)T a (0, —L, X)K~(X)r ba (X + /i)T a _1 (0, — L, -A)) 

(2.18) 

(II) SNP 

V(0,t,\,fi) = r\X) tr a (A+(A)r ab (A-/i)T a (0,-L,A)A7(A)T*«(0,-L,-A) 
+ A a + (A)T a (0, —L, X)K-(X)T^(0, —L, -A)r*j(-A - /x)) 
V(-L,t,A,/i) = r x (A) tr a (A+(A)T a (0,-L,A)r a6 (A-/i)A-(A)T^(0,-L,-A) 
+ A+(A)T o (0, — L, A)ff-(A)ri(-A - (0, — L, -A)) 

(2.19) 

Notice that all the boundary information, incorporated in K^, appears only at the boundary 
points x = 0, —L. We shall presently see that the bulk expression has in fact no dependence 
whatsoever on the reflection K matrix, since it is canceled by the t (X) factor. Note finally 
that the expressions derived in (I2.16l) - fl2.18l) are universal, that is independent of the choice 
of model. 
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3 The ATFT: brief review 



We are now in a position to systematically construct the Lax representation for any extension 

of the An ATFT following the scheme defined above, given any solution K of the algebraic 

boundary equations defined by the An Toda classical r-matrix. Note that the associated 

boundary Hamiltonians have been extracted through the asymptotic expansion of the open 

id) 



transfer matrix in [36] for sine-Gordon and in [26] for the A 2 ATFT. We shall consider in 
the following sections two particular examples, that is the prototype model of the hierarchy, 
i.e. the sine-Gordon model, as well as the Arp case. The A^ model is indeed the first 
non-trivial example of this set that may exhibit both types of boundary conditions. It is 
worth noting that in sine-Gordon the two boundary conditions coincide due to the fact that 
the model is self dual. 

Recall first the classical r-matrix associated to the generic An affine Toda field theory 



in particular is given b 

r(A) = ^ £ e, ® e, + -±- £ e^-Ci-i)^ ® e,, (3.1) 
smh(A) smh(A) 

with (eij)ki = Sik Sji. Note that the classical r-matrix (13.11) is written in the so-called principal 
gradation as in [9j [13] (see details on the gauge transformation changing the principal to the 
homogeneous gradation in [26\). We recall the Lax pair for a generic A n 1] theory [3E]: 



V(x, t,u) = -?-d x *-H+- (u e^ H E + e~^ H - - e~^ H E. e l 
2 4V u 

U(x, t ,u) = ^U-H+^ (u e^ H E + e~^ H + - e^ H E_ e^ H ) (3.2) 
2 4V u J 



<f>-H 



2A 



$, II are conjugated n-vector fields, with components 0j, 7Tj, i G {1, . . . , n}, u = e~ is the 
multiplicative spectral parameter. To compare with the notation used in [9] we set ^ = r! ^- 
{fh denotes the mass in [9]). Note that eventually in [9] both j3, in are set equal to unit. 

We also define: 

n+l n+l 
i=l i=l 

(%i are the simple roots plus the extended (affine) root, H (n-vector) and E± a . are the algebra 
generators in the Cartan-Weyl basis, and they satisfy the Lie algebra relations: 



H, E± ai 



E ai , E. 



±otiE± 



2 

-s a t -H (3.4) 

OLJ 



3 Notice that the r-matrix employed here is in fact r^* 2 with r\% being the matrix used e.g. in [111 119] 
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Explicit expressions on the simple roots and the Cartan generators are presented below. 
Notice that the Lax pair has the following behavior: 

V{x,t, -n" 1 ) = V(x,t,u), U'^u" 1 ) = U(x,t,u) (3.5) 

where * denotes usual transposition. 

We provide below explicit expressions of the simple roots and the Cartan generators for 
An [39]. The vectors = (a] , . . . , a") are the simple roots of the Lie algebra of rank n 
normalized to unity 

tth 
i 



Oi= (o,... ,0,-y^,y^,0,... ,o), i6{l,...n} (3.6) 

The fundamental weights /i^ = (/4 , . . . , /i£) , k — 1 , . . . , n are defined as (see, e.g., [39]). 

oij ■ [ik = T^j.fc • (3.7) 
The extended (affine) root a n+ i is provided by the relation 

n+l 



J2^ = °- ( 3 - 8 ) 



i=l 



The Cartan- Weyl generators in the defining representation are: 

E ai Cj i+1 , E^ ai Gj+1 i i &n+l 1 j ^—a n ^1 n+l 

n 

Hi = ^ V){ e ii ~ e j+i j+i) , < = 1 , . . . , n. (3.9) 

4 The boundary A± case: sine-Gordon model 

Let us rewrite the Lax operator for the bulk sine Gordon modefl, 

, s TUU i/3, mU~ l ifi, iff, , „ 

U(x, t, u) = —7i(x)a 3 + e^^ 3 a 2 e--^ 3 e" 0<T3 a 2 e^ CT3 (4.2) 

4i Ai Ai 

(jj are the 2-dimensional Pauli matrices. 



To recover the generic form (|3.2[) from (|4.2[) wc consider the following identifications 

--►/?, <6 u->-u (4.1) 



Here (14. 2p wc clearly consider the sine-Gordon, however after implementing identifications (|4.1|1 wc obtain 
the sinh- Gordon model. 
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Bearing in mind the expression for T ( 12.121) it is clear that we need to consider the formal 
series expansion of T and T _1 (u _1 ). But from the following symmetry of the Lax operator: 

U(u- 1 ,0,7r)=U(-«.-0,7r) (4.3) 

we see that 

T(u~ 1 ,4 > ,7i)=T(-u,-4 > ,tt). (4.4) 

We aim at expressing the term of order u in U independently of the fields, after applying a 
suitable gauge transformation [27] . More precisely, consider the following gauge transforma- 
tion such that 

T(x,y,u) = Q(x) f(x,y,u) ^(y), 
Q(x) = diag(fii(a;), ft 2 (x)) = e^ 3 ^ 3 (4.5) 

then the gauge transformed operator U is expressed as: 

U(x, t, u) = ^f(x)a 3 + ^a 2 - ^e~^a 2 e^ (4.6) 

where we define 

f(x,t)=ir(x,t) + 4>'(x,t). (4.7) 

Let T'{u) = T(u _1 ) then we introduce the following decomposition for T, T" as \u\ — > oo 



f(x,y,u) = (I + W(x,u)) exp[Z(x,y,u)} (I + W(y, u)) -1 , 

f'(x, y, «) = (!+ ^(x, «)) exp[Z(x, y, «)] (I + W(y, u))~\ (4.8) 

where the hat simply denotes that u — > —u, <p ~^ ~ c t ) - W, W are off diagonal matrices and 
Z, Z are purely diagonal. Also 

2(«) = E — ' w w = E ^V- ( 4 - 9 ) 

fc=-l k=0 

Inserting the latter expressions (14. 9 p in (12. ip one may identify the matrices and Z^ k \ 
Indeed, from equation (12. ip we conclude that the gauge transformed operators satisfy: 

dx 

+ wi] {D) - i] (D) w + wiJ io) w - u (0) = o (4.10) 

ax 



where U^, t&(°) are the diagonal and off diagonal parts of U respectively. By solving the 
latter equations we may identify the matrices Z, W. It is sufficient for our purposes here 
to identify only the first couple of terms of the expansions. Indeed based on equation (14.101) 
we conclude (see also [27]): 



m 



(4.1i; 



Note that the leading contribution as iu — > 00 comes from the e Z22 term. We assumed here 
for simplicity, but without losing generality, Schwartz boundary conditions at the end point 
x = —L, that is ir(—L) = <ft(—L) = and K~ oc I. We may rewrite the expression for the 
boundary operator V (let r a f, = Haf^ as: 



J 22 



22 



V(x, t, «, v) = r x {u)e z ™- z ™ [{I + W(0)) _1 fi(0)lir + (u)fi(0)(l + W(0)) 
1 + VT(a;))" 1 fi- 1 (a;)r(ni;- 1 )fi(x)(l + W{x)) 
1 + W / (x))- 1 fi(x)f(mOfi" 1 (a;)(l + W{x)) 



+ 



22 



but it is easy to show for the transfer matrix (12.131) as \u\ — > 00: 



t(u) 



3 Z22— Z22 



^(0))- 1 fi(0)i^ + (M)fi(0)(l + W{0)) 



22 



(4.13) 



(4.14) 



and finally 



V(x, t, u, v) 



;i + iy(x))- 1 fi" 1 (a;)r(Mt;- 1 )fi(x)(l + W(x)) 
[1 + I^(x))~ 1 fi(x)f(Hfi- 1 (x)(l + W{x)) 



J 22 



22 



(4.15) 



Again using the ansatz for the monodromy matrix we obtain from (1 2 . 1 8 [) for the end 
point x = 0: 

V(0, t, u, v) = [(1 + W(0))~ 1 fi(0)K + (u)fi(0)(l + W(0)) 
1 + l^(0))" 1 fi(0)K + (u)r(mr 1 )fi(0)(l + W(0)) 
1 + #(0))~ 1 fi(0)f(utO^ + (u)fi(0)(l + W(0)) 

J 22 



22 



(4.16) 



The r-matrix is given in (13.11) and we consider below two cases with non-diagonal and 
diagonal i^-matrix respectively. 



3 Note that in this particular case the SP and SNP boundary conditions coincide because: 



r 12 (A) = V x r^(-A) V 1} V = antidiag(l, 1) 



(4.12) 
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4.1 Non-diagonal A'-matrix 



We shall first examine the case with the generic non-diagonal i^-matrix 

A' + (A) = - sinh(A + z£)e u + - sinh(-A + i£)e 22 + x + sinh(2A)ei 2 + x~ sinh(2A)e 2 i (4.17) 

K K 

£, x ± are a priori free independent boundary parameters. The next step is to expand 
expressions (14.151) . (14.161) in powers of W 1 , and identify the first order term of the expansion. 
Taking into account the expansion of W as well as (\u\ — > oo) 

K+(u) ~ A +(0) + u~ l K +{l) + 0(u~ 2 ), riuv- 1 ) ~ r (0) + u~ l r w + 0(m" 2 ), 
r(uv) ~ r (0) + u~ l f {1) + 0(u~ 2 ) (4.18) 



where we define: 

K+<® = x + e 12 + x~e 21 , K+W = — e n - 



2 

r (0) = ^ eij <g> eu, r (1) = 2t> (e 12 (g> e 2 i + e 2 i <8> e 12 ), r (1) = 2tT 1 (e 12 <g> e 2 i + e 2 i <8> ei 2 ) 
i=i 

(4.19) 

we may expand Y(u, v) in powers of u~ l . Multiplying the resulting expression by a factor *p 
we obtain at first order: 

V(x,t,v) = —<j)'(x,t)(T 3 + fi(x,t) a 2 n~ l (x,t) + fi _1 (x,t) ct 2 Q(x,t) 

Ai Ai Ai 

(4.20) 

We see that the operator Y(x, t, v) at any point i^O coincides with the bulk operator (13. 2)) , 
consistently with the fact that coincides with the bulk sine-Gordon boundary Hamiltonian 
[4"0j I3"6] except for x = 0: 



o 

H = I dx 



2 



+ ^„WS>.^»<!!) ( 4.2i; 



Recall that in [36] the latter Hamiltonian was obtained as the first order term from the ex- 
pansion of the generating function t(u) as \u\ — > oo, assuming Schwartz boundary conditions 
at x = —L. The boundary parameters P, Q are related to the parameters £, k of the K 
matrix as: 

, _ , 'i , '£ _|_ , 'i 

P = — , Q = ■ (4.22) 

We should stress that the constraint x + = —x~ here arises by requiring that the expansions 
of the transfer matrix as iu — > oo and iu — > — oo provide the same Hamiltonians (again, 
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an analyticity condition at infinity). Such a requirement leads also to the cancelation of 
boundary terms proportional to 4>'(0). 

It is clear that the bulk V-operator is independent of the choice of i^-matrix. Expanding 
carefully the boundary expression (12.191) and multiplying the result with a factor g we obtain 
at the boundary point: 

V (6) (0, t, v) = V(0, t, v) + AV(0, t, v) , where 

AV(0,M) = -^(0)a 8 - ?(£e¥«°> + ^e-^)) ff3 (4.23) 
At 8 V2k 2k / 

V(0, t, v) is provided by the bulk expression (I4.20p . Note that all the boundary information 
is incorporated at the boundary point x = 0. The equations of motion and the corresponding 
boundary conditions emerge in this Lax formulation from the zero curvature condition. The 
zero curvature condition for the 'bulk' Lax pair yields the familiar equations of motion for 
the sine-Gordon model. We should note that analyticity requirements on the boundary Lax 
pair leads to extra constraints among the boundary parameters, i.e. x + = —x~ (see also 
|36j). The boundary operator found here is associated to the boundary Hamiltonian (14.211) . 
Note that a particular choice of diagonal K + matrix leads to discrepancies between the two 
descriptions (Hamiltonian vs Lax pair). This suggests that one has to consider as a starting 
point a generic solution of the reflection equation with several boundary parameters, which 
may then satisfy further constraints dictated by certain consistency requirements. 

The relevant boundary conditions are obtained by considering the zero curvature condi- 
tion (12.31) at the point x = 0: 

i](0,t,v) - -^-Y^(x,t,v)\ x=0 + [V(0,t,v), V®(0,t,vj\ = 0^ (4.24) 



U(0,M)-li m ^M^M^M + U(0 , M ), V(0,M) + AV(0,M) 



■ 0. 
(4.25) 

Explicit expression of the derivative of V at x = in (I4.25P indicates 

AV(0,t,u) = 0, (4.26) 

in order to eliminate a potential uncompensated divergence due to AV. 

Finally from the 'bulk' zero curvature condition and from the later expression the fol- 
lowing equations of motion and mixed boundary conditions are entailed: 

2 

171 

<f)(x, t) - <f)"{x, t) = -— wx(J3<f>(x, t)) 

M0) = |^cos(e+^(0)), (4.27) 
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which of course coincide with the equations of motion found in [10, [36] (recall also the 
identification of boundary parameters (14.221) ). In a consistent way the boundary conditions 
are obtained exactly from the Hamiltonian through 

= {H, 0(x, t)}, = {H, *(x, t)}, xe i-L, 0] (4.28) 

by noticing that the contribution containing the term 0(0) in the Hamiltonian yields a 5(0) 
term in the equations of motion for ir(x) since {(j)(x), tt(v)} = S(x — y). Elimination of this 
term yields exactly the boundary conditions (14.271) . 

Requiring cancelation of the AV(0) term is equivalent to requiring that the formal series 
expansion in u~ l coincides for V(0, t, u, v) and Y(x,t,u,v) at x — ► 0. Indeed the technical 
origin of AV(i = 0) is the non-commutation of limits x — > and u — > oo, in particular in 
e z . If these limits are required to commute then V(x,t,u) has its analytic behavior in x, u 
continued to the limit x — 0, which may suggest that it can be analytically continued "be- 
yond" the boundary. This may in turn be a relevant consistency condition in implementing 
the notion of "gluing" different boundary systems. We have thus established a straightfor- 
ward and elegant way to extract the associated boundary conditions from the zero curvature 
condition. 



4.2 Diagonal A'-matrix 

We shall now consider the diagonal .fT-matrix [4T)l W\\ 

K + (X) = sinh(A + i£)e n + sinh(-A + i£)e 22 . (4.29) 

In particular we shall be mostly interested in the degenerate limit where z£ — > oo. 

The relevant boundary Hamiltonian, obtained from the first order term of the expansion 
of lnt(A), is given by: 



o 

H — I dx 



1 . o ,/os rn 



2 



2 ^costf+f^O)) 



[2 y + V (! - cos»)j + /(0) sin - + |— . (4.30) 



The boundary contribution in ( 14.30[) is not identical with the k = limit of the boundary 
conditions in (14.211) . One needs to normalize the i^-matrix as uK(n — > 0) to get a consistent 
^-expansion, hence 7i in (I4.30P picks boundary contribution from higher orders in (I4.18p . 
By requiring the boundary term, proportional to 0'(O) to disappear we obtain the following 
constraint 

cos(£ + ^>(0)) = 0, (4.31) 
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which of course may be seen as the boundary condition to the associated equations of motion, 
as we shall see below. 

The next step as in the previous case is to expand expressions ( I4.15p . (14.161) in powers of 
it -1 , and identify the associated V-operator from the first order term of the expansion. As we 
have seen the bulk V(w, t>)-operator is independent of the choice of boundary conditions, i.e. 
the i^-matrix and is given at any point x 7^ by (13. 2p . Expanding carefully the boundary 
expression (12. 19ft and multiplying the result with a factor ^ we obtain at the boundary 
point: 

V (fe) (0,M) = 

8 # Vfi?(O)^(O)0'(O)<r 3 + ^-fii(0)fi 2 (0) (v(y + e 21 - y~e l2 ) + v~\y~ e 2X - y+e 12 ] 



(4.32) 

where A = y + fi?(0) + j/~n|(0), and y ± = ±e ±l «. 

By requiring V(0) = V^(0) we obtain the corresponding boundary conditions. Finally 
from the 'bulk' zero curvature condition and from the later expression the following equations 
of motion and mixed boundary conditions are entailed: 

2 

171 

0(x, t) - 0"(x, t) = -— sin(/30(x, t)) 

cos(^ + ^(0)) = 0, (4.33) 

which of course coincide with the equations of motion found in [10, [36] (recall also the 
identification of boundary parameters (14. 22ft and the boundary conditions found earlier. 
Notice that the obtained boundary conditions are easily obtained from the generic situation 
described in (I4.27P by simply setting the non diagonal contributions to zero. 

We are however mostly interested in the case where the .fT-matrix is degenerate (in the 
homogeneous gradation). Consider for instance the situation where K(X) = diag(e A , e~ A ). 
The Hamiltonian and boundary V-operator in this case are given respectively by: 




dx 



1, 2 , 2 , m 



2 



7T 



2 



+ — (1- cos/30) +-0'(O) (4.34) 



and 



Y ib) {0,v) = jQ 1 (0)Q2 1 (0)^e 21 -v^e^. (4.35) 

The boundary conditions emerging from the Hamiltonian are: 4>'(0) = whereas requiring 
V0>)(0,v) = ¥{0,v) in addition to the field space derivative being zero one more constraint 
is obtained: 



(4.36) 
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which of course is also automatically obtained from the boundary conditions found previously 
in the full diagonal case at = 0. Note that there is no way to trace the extra constraint 
(I4.36P from Hamiltonian point of view although we have to note that such a constraint is 
not incompatible with the Hamiltonian. 

To conclude we note that in the degenerate case some important information is automat- 
ically lost when considering the Hamiltonian description. More precisely, in the degenerate 
case there is no £ dependence so constraints of the type (14.361) disappear when examining the 
boundary conditions from the Hamiltonian viewpoint. Whenever we pass from the most gen- 
eral situation to some special situation some information is lost and inconsistencies between 
the two descriptions arise. This of course happens only when the i^-matrix possesses several 
boundary parameter and some of them are set to zero or to infinity. We shall examine in 
the following section a similar situation for the next model of the hierarchy, the A 2 theory, 
and we shall see that the arising inconsistencies may be explained in the same spirit. 



5 The boundary A^f 1 case 



We come now to the second member of the hierarchy and the first model of this class 
exhibiting both types of distinct boundary conditions SP and SNP, i.e. the A 2 model. In 
this case we have: 

ai = (l, 0), a 2 = (--, — ), a 3 = {--, — — ) (5.1) 
define also the following 3x3 generators 

E 1 = Eti = e X 2, E 2 = Et 2 = e 2 3, E 3 = Et 3 = -e 31 . (5.2) 
The diagonal Cartan generators H\ 2 are then: 

H\ = i(en - e 22 ), H 2 = ^=(en + e 22 - 2e 33 ) (5.3) 

Let T'(x,y,u) = T(a;,?/,u _1 ) and V(x,u) = U(x,u _1 ). Following the logic described 
previously (see also [27]) for the sine-Gordon model, we aim at expressing the part associated 
to E + , E- in U, U' respectively independently of the fields, thus we consider the following 
gauge transformation: 

T(x,y,u) = f(x,y,u) ^(y), 

T'(x,y,u) = n~\x) f'(x,y,u) Q(y) (5.4) 
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where we define 

Q{x) = diag(fii(z), n 2 (x), fi 3 (aO) = e ^ x) - H . (5.5) 
From equation (12.1 \ the gauge transformed operators U, U' are expressed as: 

U(z, t, u) = Q'\x) U(x, t, u) n(x) - Q-\x) 



dx 

U'(x, t, u) = Q(x) U'(x, t, u) n-\x) - Q(x) (5.6) 

dx 

After implementing the gauge transformations U, U' take the following simple forms: 

U(x, t,u) = ^S-H + ^ (uE + + -X_) , U'(x, t,u) = ^-H+^ (uE_ + -X+) (5.7) 
2 4 V u / 2 4 \ u / 

where we define: 

$ = n - d x $, 3 = n + <9,.$, X_ = e~^ H E_ e^ H , X + = e^ H E + e~^ H (5.8) 

T, U also satisfy (I2.ip . and J, J are vectors with two components fj, fj, f G {1, 2} 
respectively. 

Consider again the ansatz (j4.8p for T, T' as \u\ — > oo. As in the previous section inserting 
expressions (14. 9 p in (12.11) one then identifies the coefficients Wjj and Z^' . Indeed from (12.11) 
we obtain the following fundamental relations: 

-=#) + (P W)W 
dx 

+ wijM _ fjW W + W (i](o) w) (D) _ 0(O) _ (fj(owo) = g (5 9) 

ax 

where the superscripts O, D denote off-diagonal and diagonal part respectively. Similar 
relations may be obtained for Z, W, in this case U — > U'. We omit writing these equations 
here for brevity. 

It will be useful in what follows to introduce some compact notation: 

^ ■ H = diag(a, 6, c), ^-H = diag(a, 6, c), e^* = Ti . (5.10) 

Explicit expressions of a, 6, c and 7 j are given by: 

P,fi h . , 0, fi f 2 , £ fa 

7l = e^, 72 = e ^-^+^), 73 = e /»(-^i-^). (5.11) 
apparently a, b, c are defined in the same way as a, b, c but with fj — > fj. 
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The computation of W, W is essential for what follows. First it is important to discuss 
the leading contribution of the above quantities as \u\ — > oo. To achieve this we shall need 
the explicit form of Z^~^: 



Z^\x,y) 



m(x — y) 



/ vn 
e 3 



e 3 



1 7 



m(x — y) 



e 3 



\ 



e 3 



1 7 



(5. 



From the formulas (15.91) the matrices W^ k \ W^ k \ Z^ k \ Z^ may be determined. In partic- 
ular, we write below explicit expressions of these matrices for the first orders, which will be 
necessary in the subsequent sections (see also [26J): 



W {0) = W {0) 



( ef 1 \ 

ef -1 

7 

\ 



2iir i/K_ 

e 3 e 3 



' e 3 a c 
-a b 
\ ef c -b j 



4 



/ -b -a 
—e~~zb — c 
y a — e~3~c y 



(5.13) 



For computing the boundary conserved quantities, energy and momentum, we shall in ad- 
dition need the following expressions: 



dx 



e 3 m i \ 4e s 4e 3 

(7i + 72 + 7s) + -= (« - c') + — (a 2 + b 2 + c 2 ^ 



(i) 

22 



3 4 

e s m 



3m 



6m 



3-4 (7l + 72 + 73) + l^ (6, - a ' ) + ^ 



2 1 6 2 + c 2 ] 



(i) 
ii 



dx 



e 3 m . 

-3" "J (7l +72 + 73j 



4e 3 /W -A 4e 3 /-2 ?2 -2\ 

3m 6m 



rfZo2 e 3 m 4 e 3 f 4e 3 . 

-(7i + 72 + 73) + -5 (6' - c') + 



3 4 



3m 



6m 



a 2 + 6 2 + c 2 ). (5.14) 



5.1 SNP boundary conditions 

We shall focus in this section on the analysis of the SNP integrable boundary conditions 
in A*jp ATFT. Comparison with some already known results [9] will validate our approach, 
which then presents the advantage of being systematically implementable once a non dynam- 
ical "boundary" matrix is chosen. The boundary V-operator in this case is given by (I2.17p . 
f)2.19p . We assume here for simplicity, but without losing generality, Schwartz boundary 
conditions at x — —L and K~ oc I (see also |26j). Taking into account the ansatz for the 
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monodromy matrix (14. 8 j) as well as bearing in mind that as u — > oo the main contribution 
for the diagonal terms comes from e Zs3 , e Zs3 (see also [26]), we conclude 



V(x, t, u, v) 



'I + u))" 1 fi- 1 (a;)r(^~ 1 )fi(x)(l + 



J 33 



+ 



(i + w t (x,u))n- 1 (xy i ( u - 1 v- 1 )n(x)((i + w(x))' l y 



33 



(5.15) 



We recall that the superscript 4l denotes transposition in the first space. Also, in the ex- 
pressions with 'hat' we simply consider $ — > — $. For further technical details we refer the 
interested reader to [26]. Note that in this case the limit u — > oo is easier to consider due to 
the expressions (15.121) . In any case, although technically more involved, one can show that 
the u — > — oo limit provides the same conserved quantities and Lax pairs. In the following 
we shall expand expression (15.151) . so we need expansions of all the involved quantities: 



riuv 1 ) ~ 



0(u 



u -if (i) + 



(5.16) 



where we define: 



.(0) 



5(1) 



£ 

i=i 



2f (e 2 i <8> e i2 + e 32 <g> e 23 + e i3 ® e 31 ) 



2v 1 (e 2 i <8> e 2i + e 32 (g) e 32 + e 13 <g> e i3 ) 



(5.17) 



From the first order of the expansion of the V-operator and after multiplying with a factor 
of — we obtain: 



V(x, t, u) = -|$'(a;, t) ■ # + ^^(x, t)E + Q~ 1 (x, t) - IT 1 ^, t)E_Q(x, t) (5.18) 

which as anticipated coincides with the bulk V-operator of (13. 2p . 

In order to obtain the explicit form of the boundary V-operator we should also review 
known results on the solutions of the reflection equation for SNP boundary conditions. The 
generic solution for the A n case in the principal gradation are given by [T3l ITS] : 

n+1 



K(X) = (ge x + ge~ x ) ]T e u + ^ f ijt 



=> n + 1 v J 1 f 



+£/< 



-A— %(t-j) 



2 = 9 



g = ±q 



i=l i>j 

1 71 + 1 77 + 1 

/y = ±9 4 



i<j 



2 4 



„ H+I 

fji = q * , K J- 



(5.19) 



In order to effectively compare with the results of [13] as well as being compatible with 
[9] we always express both r and K matrices in the principal gradation (see also |26j). 
The parameter q in the solution (15. 19[) is the parameter of the underlying quantum algebra 
(quantum case) U q (sls). It is clear that we are dealing here with the classical limit q — > 1 of 
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the later solution (15.191) compatible with the classical quadratic algebra (12.91) . In this limit: 
1, 9 ->• ±1, fa -> ±1, fa -> 1, i < j- 

We come back now to the case. Recall that if + (u) = X*(it _1 ) then A + is a 3 x 3 
matrix written explicitly as: 



F = f\2 e.21 + ^23 632 + /31 ei3, 

F = /21 ei2 + /32 e23 + /13 e3i 



where 



(5.20) 



and the coefficients g, g, are given in (15.191) with n = 2 and q — > 1. Note, as in the 
previous case for the sine-Gordon model, that the boundary case for x = has to be treated 
separately. Indeed, in this case the operator V(0, u) takes the form: 



V(x = 0) 



33 



+ 



1 + # t (0))fi- 1 (0)F + (M)fi(l + W(0)) 
(1 + I^ + (0))fi" 1 (0)A' + (u)r(mr 1 )fi(0)(l + W(0)) 
1 + iy*(0))fi" 1 (0)r* 1 (n" 1 t;" 1 )A' + (n)fi(0)(l + W(0)) 



33 



33 



(5.21) 



again using the standard procedure we obtain from the first order of the expansion and after 
multiplying with an overall factor — 

V (6) (0,t,ix) = V(0,i,u) + AV(0,t,u), where 

AV(0,t,«) = ^$'(0)-iJ+^((fi 1 (0)fi2 1 (0)/ 1 2 + fir 1 (0)fi3(0)/3i)e 11 

+ (fi 2 (0)fi 3 " 1 (0)/ 2 3 - ni(0)fij 1 (0)/ 12 )e 2 2 
- (n 2 (0)fi3 ^O)/^ + fir 1 (0)fi 3 (0)/3i)e 3 3 



(5.22) 

V is the bulk V-operator given in (15.181) . The equations of motion are again obtained from 
the 'bulk' zero curvature condition: 

,2 3 



- $(x,t) + $"(x,t) 



2^ 



£ 

i=i 



/3c*j-*(x,t) 



(5.23) 



Using the same argument as in the previous section we conclude: 
AV(0) = 



<t>'M 
0' 2 (O) 



m 



(2/ 12 ef--*(°) - foeS**-*® + fsie 



f<*3'*(0) 



faeWm + / 3ie f <*-*(°) 



(5.24) 



The latter boundary conditions agree as expected with those analyzed in [9]. 

The quantities found above are associated to the first non-trivial boundary integral of 



motion: 

-0 2 2 3 



i=i ^ i=i ^ 



Q!3*(0)j 

(5.25) 



It is quite easy to check the consistency of the Lax formulation procedure in [3], i.e. that the 
Hamiltonian above leads exactly to the same equations of motion and boundary conditions as 
the zero curvature conditions. More precisely, the equations of motion from the Hamiltonian 
are obtained through: 

d(f)i(x,t) dTTi{x,t) 

— = {H, (pi{X,t)}, — = {H, TTi{X,t)} 

x e [-L, 0] (5.26) 

The boundary Lax pair in the SNP case was also constructed in [SJ. The boundary 
conditions and the Hamiltonian extracted in [26J coincide with the ones of [9j; comparison 
between our Lax pair and that derived in [9J shows that the relevant Lax pairs are also 
the same. Although we have to note that in [9] everything is expressed in terms of 9 
and S functions given that the Lax pair is derived taking into account some overlapping 
boundary regime, which in the present framework is unnecessary. Moreover in [9J the Lax 
pair construction requires a priori knowledge of the Lagrangian, given that the boundary Lax 
pair explicitly contains the Lagrangian boundary contribution. In our formulation on the 
other hand the Lax pair is derived from first principles from the expansion of (12.1 7p . (I2.19p . 
and requires no a priori knowledge of the boundary terms in the associated Hamiltonian, 
although the relevant Hamiltonian may be simultaneously obtained from the expansion of 
lnt(A). 

This ends our discussion on the sine- Gordon and SNP ATFT. We have been able to 
systematically derive consistent explicit expressions from first principle derivation for their 
Lax equations by evaluating exactly the boundary contribution. Comparison with the case- 
by-case derivation in [9] shows the same results. Our method is therefore validated and can 
be now extended to a much broader set of models beyond the ATFT. The next example will 
be the yet untreated case of A ( 2 ] ATFT under the "soliton -preserving" reflection algebra (see 
also on SP boundary conditions within ATFT). 
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5.2 SP boundary conditions 

We shall now examine the SP boundary conditions (see |26j). We shall consider below two 
cases with the A-matrix being non-diagonal and diagonal respectively. 

5.2.1 Non-diagonal A'-matrix 

We choose for A the non-diagonal solution of the reflection equation found in [HI and 
demonstrate how this particular choice of A'-matrix contributes to the integrals of motion, 
and the relevant boundary Lax pair. We consider for simplicity Schwartz boundary condi- 
tions at the end point x = —L, whereas the right boundary is described by this A-matrix 

gang 

K+(u) = V K{u) V with V = diag(l, e% , e^) 

K{u) = (e 4A - l)g 2 + (e 4A - l)(£e 2A - a - f3)g - (£ 2 e 4A - (a + /3)£e 2A + a/3)I(5.27) 

where we define the 3x3 matrix (HJ US] : 

g = (a + (3)e n + x + ei2 + x~e 2 i, x + x~ = -a(3 (5.28) 

a, (3, £ are free independent boundary parameters, and x^ 1 satisfy ( 15.281) . To proceed with 
the expansion of lnt(A) and V in powers of u~ l we shall need the expansions of r, f, A + , 
(as \u\ — > oo): 

r(uv) ~ r (0) +u~ 1 f {1) + 0(u~ l ) 
K + (u) ~ A (0) + u~ 1 A (1) + £>(u~ 2 ). (5.29) 

where is defined in ( 15. 17ft and 

f (1) = 2tr 1 (ei2 <g> e 2 i + e 23 <8> e 32 + e 3 i <g> e 13 ), 

A(°) = -e 2 e 33 + £x + e 12 + £x-e 21 , A« = £(a + 0)e n (5.30) 

We shall consider henceforth for simplicity x + = x~ , and also set y = x^^. The integrals of 
motion follow from the asymptotic expansion of the logarithm of the open transfer matrix 
as iu — > ±oo, i.e. 

I 7 (n) T 

i nt{iu _ _oo) = y in = y Z22 ~ Z22 + y ^ (5.31) 
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More technical details on such expansions will be found by the interested reader in [26]. The 
integrals of motion obtained from the first order of the expansion are given by: 

ii = -^—(v^ + iVsn^), 

12m 

Ii = - — (V (b) -tV3H {b) ). (5.32) 
12m 

The momentum and energy are directly obtained from the above conserved quantities and 
defined as: 

= f d X j2 (** # - < <t>) + X>(o) mo) - + mi(o) 

J L i=l i=l P 

'yfti(0)ft 2 (0) ( - 4a(0) + 6(0) - 4a(0) + 6(0)) - £ 2 O*(0) (c(0) - 6(0) + c(0) - 6(0) 
yfii(0)fi 2 (0) f c(0) - a(0) - 6(0) + a(0)) + (c(0) - 6(0) - c(0) + 6(0)' 



/3 2 A 



/3 2 A 



(5. 



where A = yili (0)^(0) — ^ 2 O|(0), and f2 and a, 6, c, a, 6, c have been previously defined. 
As already pointed out in [26] the two boundary cases exhibit essential differences: in SNP 
the c-number if-matrix contains no free parameters, and consequently no such parameters 
occur in the deduced integrals of motion. In the SP case on the other hand the i^-matrix 
contains free parameters, which explicitly appear in the boundary integrals of motion. More- 
over, in the SNP case, only the boundary Hamiltonian belongs to the family of commuting 
quantities, whereas now both boundary Hamiltonian and momentum turn out to be con- 
served quantities. Note that we are here in a situation where the terminology "boundary 
effects" is misleading, since it suggests that no conserved momentum could exist due to the 
presence of a physical boundary. As we have shown here this is not true, indicating that the 
physical interpretation of this field theory is more subtle and has to be further explored. 

We shall now derive the associated boundary V-operator. From the two different limits 
(iu — > ±oo) we obtain essentially two quantities from the first order of each expansion: 

Vf\x,t,v) = (v(x,t,v) — iVS¥(x,t,v)] from iu — > oo 

3m V / 



Y { 2 ] (x,t,v) = — [U(x,t,v) +iv / 3V(x,t,w)) from iu — >• — oo (5.34) 
3m V / 



where V, U is the Lax pair of the A^p theory defined in (13. 2p . The bulk operator as already 
mentioned is independent of the choice of i^-matrix. 
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The expressions of the V-operator for each end point are given below: 



V(0,t,u,v) 



.1,1 



x 



+ 



1 + W A (0))- 1 fi(0)A^+(M)fi(0)(l + W(0)) 
1 + W(0))" 1 fi(0)K + (u)r(ui;~ 1 )fi(0)(l + W(0)) 

i + t^(o))-Hi(o)f(H^ + («Mo)(i + w(o)) 



33 



33 



J e {1, 2}. 

(5.35) 



At the boundary point we get: 



- uW(o,«; 

m v 



— > OO 



-oo. 



(5.36) 



Let us focus on (I5.36p . which is associated to the Hamiltonian of the system (15.331) 
as will become transparent in the following. Note that similarly plays the role of the 
V-operator associated to the momentum of the system. We then define: 



V (6) (0, 



= ^ 2 ^j(0)yfii(0)ft 2 (0) ( c (o) - 5(0)) (e n + e 22 - 2e 33 ) 
+ 2A ( - ^ 2 (0)fi 3 (0)e 23 - yfi 2 (0)fi 3 (0)e 31 j 
- ^ 2 (0)fi 3 (0)e 32 - yfi 2 (0)fi 3 (0)e 13 N 



2A 



U {b) (0,i;) 



+ 
+ 

+ 



A 



^£ 2 fi*(0)yfli(0)fi 2 (0) (a(0) - 6(0) + a(0) - 6(0)) (e n + e 22 - 2e : 

a(0) + 0(0)) fyfii(0)Q 2 (0)(e la - e 22 ) - 2£ 2 fi 2 (0)(e n - e 33 



2A 



mv 
2A 



2A 



2yfi 2 (0)e 12 - £ 2 fi 2 (0)fi 3 (0)e 23 + yfi 2 (0)fi 3 (0)e 31 



2yfi 2 (0)e 21 - £ 2 fi 2 (0)fi 3 (0)e 32 + yfi 2 (0)fi 3 (0)e 13 



We shall focus on the equations of motion emerging from the Hamiltonian via (15.261) . 
We expect that identical equations arise from the zero curvature condition. Indeed the bulk 
part gives rise to the equations of motion that coincide with the familiar ones (see equations 

dE2D). 

From the Hamiltonian derivation the boundary conditions arise by requiring the boundary 
terms (proportional to <^(0)) to vanish, yielding: 



$'(0) = yfii(0)Q 2 (0) = -rfia(O). 



(5.37) 



From the zero curvature condition however, as analyzed in section 2, we require V(0) = 
V^(0), and we end up with an extra constraint in addition to the ones (15.371) emerging from 
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the Hamiltonian, i.e. 



fii(0)fij x (0) = 0. (5.38) 

Although this extra constraint is compatible with the Hamiltonian description it is neverthe- 
less missing when analyzing the boundary conditions from the Hamiltonian point of view. 
This 'missing' information may be associated to the fact that the i^-matrix we have chosen 
( I5.27P is not the most general matrix to start with, and exhibit some degenerate or at least 
non-generic behavior, as we shall comment in the final section. Finally one may derive the 
equations of motion from the boundary V-operator associated to the momentum V along 
the lines described above. 

To conclude we have here also been able to explicitly derive exact expressions of the 
boundary Hamiltonian and momentum -both conserved quantities for the particular bound- 
ary conditions- as well as the associated boundary Lax pairs. 



5.2.2 Diagonal i^-matrix 

We shall now focus on diagonal degenerate solutions of the reflection equation [JT] given by 
the following expressions (in the principal gradation): 

l n 

K m (X, f) =sinh(A + iOe" A ^e~^ (j " 1) e ji + sinh(-A + iOe A ^ e"^ '"" 1 ^ (5.39) 

3=1 j=l+l 

2A 

(recall u = e n + x ). As we have seen from the previous section on the sine-Gordon degenerate 
iT-matrices give rise to inconsistencies. Some information is usually lost from the Hamil- 
tonian point view when taking the limit from the generic non-degenerate K matrix to the 
degenerate one, which is a special case of the generic solution. To obtain the i^-matrix in 
the homogeneous gradation one implements a simple gauge transformation (see e.g. [26]). 

In the A2 case we end up with two types of diagonal boundary matrices corresponding 
to the two possible values 1 = 1, 2. We shall consider an example here to demonstrate how 
the particular choice of boundary i^-matrix contributes to the integrals of motion. It is 
practical for the following to consider a non-trivial left boundary described by Km, and a 
right boundary described by the i^( 2 )-matrix, i.e. 

K + (u, t) = Km{u-\ £+), K~(u, f") = K m {u, £")• (5.40) 

We now proceed with the expansion of V in powers of u~ l . The bulk part, for 0, —L, 
is given in (15.341) . 
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To obtain the generic results with the least effort it is practical to consider the two 
different types of i^-matrices to each end of the theory. The expressions of the V-operator 
for each end point x b = 0, — L are given by 

i 



V? ) (x b , t,v) = - (V^ (x b , t, v) - iVsV® ( Xh t, v)) 
m V / 



for iu — > oo 



Y (b \x b ,t,v) = —(i] {b) (x b ,t,v)+iV3V ib) (x b ,t,v)^j for iu — > — oo. (5.41) 

We are mostly interested in V^ , which we shall associate to the Hamiltonian of the system 
as will become transparent in the following. Note that plays the role of the V-operator 
associated to the momentum of the system. We then have 

VW(0,v) = j^nf^tve^-v-'e^ 

Y( h \-L,v) = ^(-L^i-L^ve^ + v- 1 ^. (5.42) 

Just for the record we also give: 

UW(0,t;) = ^ e -^ + ^(0)fi3 2 (0)(e 3 3-en) +^(0)fi 3 " 2 (0) (^3 + ^32) 

U {b) (-L,v) = ^e-^ + fi 2 (-L)fi3 2 (-L)(e 33 -e 11 ) + ^ 2 (-L)fi 2 ~ 2 (-L) (t;e 12 + rt 21 ) . 

The integrals of motion emerging from the first order of the asymptotics of the transfer 
matrix as iu — > ±00 are given by: 

Ii = = -^-(V^+iVSH {b) ), 
12m 

Ii = -J?—(V {b) -iVSH®). (5.43) 
12m 

The momentum and energy are directly obtained from the above conserved quantities and 
defined as: 

7>W = f dxJ2 Ui # ~ < + X>i(0) MO) + |a a ■ H(0) + l^ e - 2i « + fi 2 (0)fi3 2 (0) 
i=i i=i ^ ^ 

2 

- £ ^(-L) fc(-L) - ■ H(-L) + ^e- 2 ^fi 2 (-L)^3 2 (-L) 

^ = J\x{j^tf + tf) + ^j^e^ (5.44) 

Naturally the two boundary cases are distinguished; in SNP the c-number i^-matrix contains 
no free parameters, and consequently no free parameters occur in the entailed integrals of 
motion. In the SP case on the other hand the iT-matrix contains free parameters, which 
explicitly appear in the boundary integrals of motion. 
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As mentioned earlier upon the 'boundary' quantity is associated to the Hamiltonian 
HS h } . We shall focus on the equations of motion emerging from the Hamiltonian via (15.261) . 
The same equations should emerge from the zero curvature condition. Indeed the bulk part 
gives rise to the equations of motion that coincide with the familiar ones. Note however that 
from the Hamiltonian the boundary condition are: 

$'(0) = = (5.45) 

whereas from the condition V^(xf,) = Y(xb) in addition to (15.451) extra constraints involving 
the boundary field appear, more precisely: 

fi^ 2 (0) n 2 3 (0) -> ^i 2 {-L) ft 2 2 {-L) -> 0. (5.46) 

Again as in the non-diagonal case a discrepancy between the two descriptions is apparent 
due to the degenerate nature of the i^-matrix. 



6 Comments 

We have analyzed via the boundary Lax pair formulation two distinct types of boundary 
conditions. The SP case presents a particular interest given that certain subtle technical 
points arise requiring further clarification. Let us briefly comment on the difficulties emerging 
in this context. 

We see that when the i^-matrix possesses several free boundary parameters, discrepancies 
between the Hamiltonian and Lax description, which however are expected to be equivalent, 
are seen to emerge when some parameters are equal to zero (e.g. for SP boundary condi- 
tions in case and even in sine-Gordon). This is an intriguing point and its resolution is 
probably associated to defining a notion of "appropriate choice" of the c-number representa- 
tion (if-matrix) of the reflection equation. A generic non-diagonal solution of the reflection 
equation is required as a starting point. For instance in sine-Gordon model, amongst all 
non-diagonal solutions (I4.17P the only "legitimate" solution is the one with x + = —x~ . Such 
a restriction is here dictated by the requirement of a consistent asymptotic behavior of the 
generating function. 

It is worth noting that diagonal ^-matrices we have considered here as well as non- 
diagonal ones of the type [42] in (SP) have only two distinct eigenvalues (homogeneous 
gradation) hence their spectra are doubly degenerate. Similar discrepancies occur as we have 
seen in the context of sine-Gordon model when choosing to consider K oc I (homogeneous 
gradation). Such a matrix has just one distinct eigenvalue so it is doubly degenerate, leading 
again to inconsistencies. In the example we considered here in the A^ case even though the 
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non-diagonal matrix K matrix is not degenerate (homogeneous gradation), the matrix g in 
(15.27ft has a zero eigenvalue, and presumably this is the point that creates the problem in 
this case. 

In a more algebraic framework, all solutions of the reflection equation -the solutions of 
the quantum reflection equation satisfy the classical reflection equation as well- considered 
in the present work are representations of the cyclotomic Hecke algebra (see e.g. [^3lH3lH5] ). 
The generators g\, go, I = 1, . . . , iV — 1 of the cyclotomic algebra Cjy satisfy the following 
set of constraints: 

9i gi+i gi = gi+i gi gi+i, gi go gi go = go gi go gi, 

n 

(gi-q)(gi + q~ 1 )=o, H(g -Q = o 

a=l 

[gu g m ] = 0, \l-m\> 1, [g , gi]=0, I > 1 (6.1) 

One expects that the most generic solution for the A^li case should be expressed in terms 
of representations of the generator g with all £ Q ^ and £ a ^ £p V I ^ k, i.e. (see also 

(E2ZD) 

n-l 

K(X)=^Ca(X)^ (6-2) 

a=0 

the general solution should be thus dictated by the rank of the algebra. Note in particular 
that for the sine-Gordon model n = 2 one recovers the boundary Temperley-Lieb algebra 
(see e.g. [12]). In this context all solutions that give rise to inconsistencies are special in 
the sense that are either degenerate £ Q = £p or correspond to a case with at least one zero 
eigenvalue £ Q = 0. 

To summarize: special cases of ^-matrices give rise to inconsistencies, hence one needs 
to consider the most general possible solutions of the reflection equation with distinct inde- 
pendent boundary parameters. In the SNP case no extra free boundary parameters appear 
and no extra constraints among the boundary fields occur. In the A { 2 ] SP case we conjecture 
that any generic (non-degenerate) non-diagonal solution with free boundary parameters will 
be appropriate. For the moment we have no such a generic matrix at our disposal, but 
the inconsistencies arising give us a strong hint that there should exist solutions with more 
boundary parameters. We thus conjecture that the K matrix ( 15.271) will turn out to follow 
from some yet-to-be found general solution via a limit process. This procedure is causing loss 
of information in the Hamiltonian analysis (possibly through a subtlety in the formulation of 
the exchange between this A'-matrix limit and the asymptotic expansion limit which yields 
the boundary contributions), giving rise to the observed inconsistencies. We have thus a 
strong motivation to systematically search for more general solutions in the SP case. 
We shall further pursue this significant issue in a separate publication. 
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